Geometrical approach to light in inhomogeneous media 
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Electromagnetism in an inhomogeneous dielectric medium at rest is described using the methods 
of differential geometry. In contrast to a general relativistic approach the electromagnetic fields 
are discussed in three-dimensional space only. The introduction of an appropriately chosen three- 
dimensional metric leads to a significant simplification of the description of light propagation in an 
inhomogeneous medium: light rays become geodesies of the metric and the field vectors are parallel 
transported along the rays. The new metric is connected to the usual flat space metric diag[f ,f ,f] 
via a conformal transformation leading to new, effective values of the medium parameters e and fi 
with ijl = 1. The corresponding index of refraction is thus constant and so is the effective velocity 
of light. Space becomes effectively empty but curved. All deviations from straight-line propagation 
are now due to curvature. The approach is finally used for a discussion of the Riemann-Silberstein 
vector, an alternative, complex formulation of the electromagnetic fields. 



I. INTRODUCTION 



During the past few years much research has been done in the field of analog laboratory models of general relativity. 
^ ■ A particularly fruitful approach was the discussion of the propagation of signals in moving media. In an inhomo- 
r — \ geneously moving medium, light Q or sound signals deviate from a straight path creating the impression of an 
. attractive force drawing them towards some point. This effect has led to the notion of optical and acoustical black 
holes, velocity profiles where all signals coming sufficiently close to some point are so strongly attracted that they 
, fall into it. From the theoretical point of view the most intriguing result of these investigations was the realization 
that these effects can actually be described in geometrical terms; an effective space-time metric creates a space-time 
[ curvature and thus leads to the deviation from the straight path. Although the idea of creating table-top experiments 
being models of astronomical objects triggered a particular interest in the physics of moving homogeneous media, 
they are by no means the only systems where differential geometry can be applied to describe the propagation of light 
in a dielectric. 

The system described in the present paper is an inhomogeneous dielectric at rest, i.e., a dielectric with electric 
permittivity e and magnetic permeability ii depending on position. The best known everyday example of such a 
medium is the air above a hot road which reflects light coming from the sky. Here the effect is due to the different 
temperatures in adjecent layers of air which lead to a variable density. A more extreme example is the interface 
5J] ' between air and glass where the index of refraction n = ^/Jli varies discontinuously, but this kind of system will not 
' be considered here. 

The propagation of a light ray in an inhomogeneous medium has usually been described in terms of Fermat's 
principle: light moves along paths with the shortest optical length, i.e., the paths that give the shortest transit time. 
The fact that light rays extremize some measure of length suggests that they might be geodesies of the corresponding 
metric. That this actually is the case with a metric of the form ds^ — v?{dx^ -f dy^ + dz^) has been shown by 
Bortolotti in his 1926 paper H). The propagation of the field vectors along the ray has been addressed for the first 
time by Rytov [Q who discussed this topic in terms of classical differential geometry, calculating the rotation of the 
field vectors in relation to the principal normal of the trajectory. Kline and Kay show in their book |^ that the field 
vectors are actually parallel transported along the curve in the sense of Levi-Civita parallelism. Unfortunately, the 
authors stick to a mathematical notation as used in the differential geometry of curves which makes the calculations 
rather cumbersome and conceals the physical meaning of the results. Finally, Solimeno, Crosignani, and DiPorto 
mention that it is possible to introduce a metric tensor and identify the rays as geodesies in this metric. Furthermore, 
they represent the ray equation in the form of a geodesic equation in terms of Christoffel symbols. They do not, 
however, describe the fields and the field vectors in this formalism. 
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The aim of the present paper is to give an account of how electrodynamics in an inhomogeneous medium can be 
represented in an alternative metric. We will extend the discussion beyond the ray equation discussed in |6|. With a 
few exceptions all calculations will be performed in the convenient notation of the Ricci calculus as used in the general 
theory of relativity. This keeps the calculations simple, gives them a clear meaning and makes them easily accessible 
to anyone acquainted with the methods of general relativity. We will confine our attention to metrics connected to 
the standard metric by a conformal transformation, i.e., they differ by a scalar factor. Some of the calculations are 
valid in the present form for a more general class of metrics, but as these metrics are of no physical significance for 
the systems considered here they will not be discussed. The more general system - the anisotropic medium - requires 
still more sophisticated mathematical methods. 

The paper is built up as follows: In sec. || we introduce the general conformally transformed metric and show how 



the fields /i, e and the wave equation are transformed under a conformal transformation of the metric. In sec. Ill 
the approximation of geometrical optics is introduced and the wave equation is rewritten as a series in powers of the 
inverse wave number. Imposing on the lowest order term the condition that the trajectory has to be a geodesic leads 
to the metric introduced by Bortolotti. In sec. |^ the propagation of the field vectors is discussed together with the 
intensity flow. The wave equation in the new metric is given. Finally in sec. the Riemann-Silberstein (RS) vector 
in an inhomogeneous dielectric is discussed. This concept was introduced by Riemann 0| and its main properties 
have been discussed by Silberstein Using the name Riemann-Silberstein vector, we follow Bialynicki-Birula 
The RS-field allows for an alternative, complex, formulation of electrodynamics, and is considered a candidate for the 
photon wave function M. 



II. FIELD EQUATIONS 

In the present section we consider the electromagnetic field equations in an inhomogeneous dielectric medium. As 
usual the electric and magnetic fields in this system can be described by Maxwell's equations 

V-B==0. VxE = — — , (1) 

dT) 

V-D-0, VxH-— . 

ot 

The properties of the medium are concealed in the connection between the fields, which in our case can be written in 
the form of the constitutive equations: 

D = eeoE, B = ^^qH (2) 

where we allow the electric permittivity and the magnetic permeability to depend on position. Maxwell's equations 
(|l|) and the constitutive equations (||) allow us to derive the wave equation in the medium and, consequently, to 
discuss the propagation of light. Unfortunately, in general the equation obtained in this way includes many terms 
which do not have a clear physical interpretation. In this section we show that the introduction of an alternative 
space structure - a changed metric - can simplify the equations considerably and give them a clear physical meaning. 

A metric is required, whenever one wishes to define lengths of vectors. The length of a vector a can then be written 
in the form 



|a| = y/a^a^ (3) 

where 7^j^ is the metric tensor. For the reader familiar with general relativity we should emphasize that the calculations 
are performed in a three-dimensional space and that the metric tensor of flat space is the unit matrix. In (||) we 
assume Einstein's summation convention, dropping the summation symbols and implicitly summing over repeated 
indices. In a space equipped with a nontrivial metric tensor, two different versions of the same vector, the co- and 
the contravariant one, have to be considered. They are mapped into each other by the metric tensor and its inverse. 
These tensors are said to raise and lower the vector's indices because the co- and contravariant vectors are denoted 
by lower and upper indices, respectively: 

= a'' 7^^. (4) 

In ordinary three-dimensional electrodynamics there is no need for an explicit introduction of the metric tensor which 
is simply the identity matrix, and the components of the co- and contravariant vectors are equal. In the present 
paper, however, we will allow for more general metric tensors which correspond to a more general form of Maxwell's 
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equations. When Maxwell's equations are written in coordinates, the antisymmetry of the curl operator is represented 
by the antisymmetric tensor e°'^'^ . In the case of general metrics this has to be replaced by the generalized form: 

a/37 ^ J_e"/37 (5) 

with 

7 ^ det(7^^). (6) 
The covariant metric tensor 7^,^ is the inverse of 7'"': 

lu.r' = K- (7) 

Following the procedure used by Landau and Lifshitz [ |lO[ one can rewrite Maxwell's equations in the form 

r^-^^pE, + ^ = 0, r^-^^pH, - ^ = (8) 

and 

V„D" = 0, V^S" = 0. (9) 

The symbol V/5 denotes the covariant derivative, the generalization of the partial derivative d/dx^ to curvilinear 
coordinates or curved spaces. In coordinates, the covariant derivative of a vector field a^^ acquires the form 

V„a^ = |^ + r^„;3a^. (10) 
The Christoffel symbol r'^^^ can be calculated when the spatial dependence of the metric tensor is known 



~ 2^ V dxt^ dx- ) • ^ ^ 

It accounts for the fact that not only the vector components but also the coordinate vectors themselves depend on 
position. Consequently, the covariant derivative of a scalar field is equal to the partial derivative. Details can be 
found in any book on general relativity, e.g., [Tl|Jl^ . In the present paper we will avoid the application of explicit 
coordinates, since our goal is to show the fundamental geometrical structure experienced by light propagating through 
an inhomogeneous medium. Obviously, as soon as one is interested in a concrete example, one has to choose a system 
of coordinates. 

Maxwell's equations in the form (||) and have been introduced formally as an alternative to the standard form 
(|l|). They are obviously correct when the standard metric is used, but one has to check if they remain valid in the case 
of alternative metric tensors. The metric tensors we consider here are connected to the flat metric via a conformal 
transformation, i.e. 

= ^^Sf,, (12) 

with n being some differentiable positive scalar-valued field. 

When changing to a conformally equivalent metric it may become necessary to introduce also transformed fields of 
the form 

Da = ^''Da (13) 

where s is an integer called the conformal weight. In order to calculate it one has to use the explicit form of the 
covariant derivative (|l^) with the Christoffel symbol in the form 

Eq. then acquires the form 

r^^aDp = n-^j'^^id^nwp ~ n^ripD^) (15) 
= n-^{{s - 2 + d){n''-^D'^df3n) + n^daO"^) (16) 
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where d is the dimension of the space, which in our case is 3. For eq. to be equivalent to the usual Maxwell equation 
the first term in (|T^) has to vanish, leading to the result 

(17) 

which therefore implies 

Da = n~^Da and = fl-^D" (18) 

with an analogous relation for the B-ficld. Because the Christoffel symbol is symmetric in the two lower indices we 
can rewrite eq. (||) as 

±e-^^^ + ^^0 (19) 

and, taking into account that ^ = fl'^, one gets 

= Ea and E'^ = Q-^E" (20) 
with the equivalent result for the _ff-field. These results lead to the alternative constitutive equations 

Ba = n^^Ba = n^^fl^loHa = l-lfloHa (21) 

Da = n-^Da = n-^eeoEa = n^^eeoEa, (22) 
valid in the transformed metric. This allows us to define the permeabilities of the system in the new metric as 

In order to arrive at the general form of the wave equation we calculate as usual the curl of the first equation in (|^) 
applying the covariant curl operator with rjafi-y = y/lGafi-y, and get 



This leads to the equation 



- 5;<S^) V'-V^^^ + ^ry^p^V'-B" - 0. (25) 



Here we applied the fact that the covariant derivative of the modified antisymmetric tensor ri°'P'^ vanishes. This can 
be easily seen, taking into account that the covariant derivative of the antisymmetric tensor e°'^'^ does not vanish, 
but rather leads to the result 



y a/37 ^ r-l^e"^^ = ^!^^/2e"^T (26) 
' oxP 



where in the last step the well-known result for the contracted Christoffel symbol has been applied [0. Note the 
order of indices in the antisymmetric tensor: 

Vc.paVPB'^ =. -rj^paVPB" (27) 

with the last term corresponding to the curl of B°'. Thus we get for (E5h 



V^VpE, - V^V,^^ - ^^'q.pJ/PB - 0. (28) 

In order to be able to apply the divergence equation (|^), one has to interchange the covariant derivatives in the 
second term. But, as one easily sees from the explicit form in eq. (|lO|), covariant derivatives do not usually commute. 
The commutator of covariant derivatives is connected to the curvature of space, e.g., when applied to a vector the 
commutator gives 
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[V„,V^]a^ = i?%„/3a^ (29) 



with R^vap being the Riemann curvature tensor [ |liyi2f| . Note that the conformal change of the metric tensor may 
create a curvature, and the fact that the considered space is flat in the standard metric does not imply that this 
remains so after the transformation of the metric. In our case, the commutator gives 

[V„,V^]^^ = i?a^^'' (30) 

where 

Rati = R^ aPu (31) 

is the Ricci tensor. Finally, we get the wave equation in the form: 

V^VpE^ + V,(^^V^ln£) - ^Q^E, + R^^E^ + {V^Ep - Vp^,)V''ln/i = 0. (32) 

One easily notes that when the flat metric is applied, this is the usual wave equation for light in an inhomogeneous 
medium, e.g. 0: 

us (9^E 

V^E + H-^-^ + (grad In /i) x (V x E) + grad(E • grad In e) = 0. (33) 



III. GEOMETRICAL OPTICS 



The physical content of the metric transformations becomes clear in the framework of geometrical optics, i.e., in 
the case of monochromatic hght moving in a medium with the medium parameters not changing significantly within 
one wavelength. This approximation allows us to introduce the notion of light rays propagating through the medium. 
The field can then be written in the form 

Ea = 4exp(i(fco5 - ujt)). (34) 

£a is the envelope of the wave; other fields like or I?" are defined likewise, w is the constant frequency and fco is 
defined as fco = w/c = 27r/Ao with Ao being the vacuum wave length. The product of fco and S - the "optical path" 
- is the spatial phase of the wave. Inserting ( ^4[ ) into the wave equation (^) and expanding in orders of 1/fco gives: 

0= -£„\If3S\l'^S + iie£c, (35) 

+ hv'^iaVfjS + iaV^VpS + £-yVaSV''\ii {/is) - £^VpSVnn jl 
fco ^ 

+ ^ (\7^^p£„ + V'^ln {fie)\/a£'y + V'^ln/iV/j^ - Vpf^V^ln/l + Rc,f,£^ 

fcg V 

One should not be confused by the fact that the expansion coefficient has a dimension. The crucial point here is that 
two kinds of spatial variation are present here: the fast variation of the phase and the slow variation of the envelope. 
The clear separation of these scales is not a result of the calculations but rather a condition for the geometrical 
approximation to be valid. This validity is reflected in the size of fco which has to be large enough for the three terms 



in (35) to be well-separated. 



Confining our considerations to the zeroth order contribution, we get the equation 

V'^^V^j^ jle. (36) 

In a homogeneous medium light rays follow straight lines, but in an inhomogeneous one this is usually not the case. 
In the case of a general metric the idea of the "straightest line" is represented by the concept of a geodesic - the line 
that "curves as little as possible" |ll[] and that extremizes the distance between two points. In mathematical terms 
the geodesic is most easily described as the line that "parallel transports its own tangent vector", i.e., the covariant 
derivative of the tangent vector field in the direction of the curve vanishes. Calling a tangent vector t" we thus get 
the condition 
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f'Vc^t^ = (37) 

for a curve to be a geodesic. In the case of propagating light, the tangent vector is the velocity vector V"5, the 
contravariant version of the wave vector V^iS, which is the gradient of the phase S. Consequently, for the light ray 
to be a geodesic the condition 

VSVo^V^S = (38) 

must be fulfilled. Note that 5 is a scalar and therefore its first derivative does not depend on the metric and two 
derivatives applied to it commute. Eq. (|3^) can be rewritten in the form 

(V"5V„5) = 0. (39) 

Since the expression in the brackets is a scalar, we end up with the condition that for the light trajectory to be a 
geodesic, V"5Vq,5 has to be constant. This is equivalent to the condition that e/i has to be constant. This condition 
allows us to find the "most natural" metric for a given medium. Using ( p3| ) we conclude that the light ray is a geodesic 
in the metric 

/n2 \ 

n2 , (40) 
J 

i.e., one has to apply the conformal transformation with 

= n^. (41) 

In ( ^ ) we denoted the metric that allows us to represent the trajectories as geodesies with a bar over the symbol. All 
objects corresponding to this metric will be denoted by bars over the respective symbols, whereas symbols without 
bars correspond to arbitrary metrics of the form ([T^). Obviously, the metric tensor in ( ^ ) may still be multiplied 
by a constant but our choice makes p. and s inverse. Note in particular that the product of p, and e corresponds to 
an effective index of refraction. Consequently, the metric introduced here leads to a constant index of refraction and 
thus to a system with a constant effective velocity of light. 
The quantity p. introduced here is equal to 

'ft. ,42, 

Together with the index of refraction n, m gives an alternative description of the medium properties. The interesting 
point is that m drops out in the calculations. Therefore it is possible to give a complete description of ray propagation 
in an inhomogeneous medium in terms of geodesies of a metric tensor. The quantity m is sometimes called the 
impedance M or the resistance M of the medium. This is due the the fact that it has the units of Ohms. 



IV. POLARIZATION TRANSPORT 



Since the electromagnetic field is not a scalar field, the knowledge of the trajectory of a light ray - the set of the 
points it passes through - does not contain complete information about the propagation. One still needs to find out 
how the field strength vectors are propagated along the ray. In order to address this question we go back to eq. (|3^), 
assume that eq. ( p6| ) is fulfilled and demand that even the first order contribution vanishes. This leads to the equation 

2V^£^VpS = -SaV^VpS - L,V^SV''ln {efl) + f ^ Vp5V''ln jl. (43) 

In analogy to (|3^ ) the vector is said to be parallel transported, along the ray when the condition 

V^£aV,}S ^ (44) 

is fulfilled, the left-hand side corresponding to the directional derivative of along the ray. In a more general sense 
the vector £a is still said to be parallel transported if 

V^SaVpS = a£, (45) 
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but in this case the length of the vector is not conserved during the propagation. One sees immediately that in order 
for the field vector to be parallel transported the second term on the right hand side of ( ^3|) has to vanish. This 
is equivalent to the condition that is constant. Not surprisingly, this is exactly the same condition as for the 
trajectory to be a geodesic. Consequently, in this case we get the equation 

2Wl^£^VpS = -4 V'^ + 4 Vp^V^ln fl (46) 



where we used the bar as introduced at the end of section III, The last term vanishes, if /i is constant. This corresponds 



to the condition that fi/n = y/e/fi is constant in the standard (or any other) metric, i.e., e and fi differ only by a 
constant factor. This means that jl and e are constant. In this non-physical but interesting case, we get for the 
constitutive equations 

E°' =D°' B" = H". (47) 

This means that the whole contribution of the medium iscontained in the metric structure of space, the light is propa- 
gating as if it were moving in empty but curved space. This result is not only valid in the geometrical approximation, 
but even allows us to write the wave equation in terms of a curved space as follows: 

ypypE^Y^ + R^^E^^ - = 0. (48) 



This case corresponds to a constant value of the quantity m introduced in (42). In the general case of arbitrary // 
and e we get for the transformed wave equation: 

Vf^VpK - f^^Q^^„ H- ^^V<,V^ln£ Vp^^V^lne = 0. (49) 

As £ = 1/m, this equation shows that in the discussion of the propagation of the field vectors in the general case both 
n and m have to be used. Since the geometrical description is based on the index of refraction n only, no complete 
geometric description is possible in the general case. 

With the field vectors being parallel transported along a geodesic, the angles between the field vectors and between 
the fields and the ray remain unchanged during the propagation along the ray. It is thus sufficient to find these angles 
at one point on the ray in order to know them for the whole trajectory. One can easily convince oneself that in our 
approximation all the angles are right angles. This can be seen from the geometrical approximation of Maxwell's 
equation by inserting the harmonic fields (jsj) into Maxwell's equations (||) and (|^) and keeping only the lowest order 
in 1/fco- We drop the tilde over the symbols when the general case is meant. We thus get 

•q"''^'^ {V pS)n^ + ££0^" = and 77"'^'' (V;,^)^.^ - /^Mo^" = (50) 

and 

r'"Vc.5==0 and S"V„5 = 0. (51) 

The last two equations show that the field vectors are orthogonal to the ray. In order to show that they are orthogonal 
to each other one has to multiply the first equation of (^) by Tia which leads to 



££of"H„ = (52) 

where the antisymmetry of the tensor rj'^P^ has to be taken into account. Finally, contracting the first equation of 
( pO| ) with £a and the second one with Tia we arrive at 

££o£"£a = fJ-^^oTi^Ha- (53) 

Note that the terms in eq. ( p3[ ) are - up to a constant - the electric and magnetic energy densities, respectively, 
showing that both densities are equal. The complete energy density of either field can be written in the form 

We = i^I?" = (54) 
= ^HcB^' = ^Ho^n". (55) 
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The transformation of the fields to the new metric as described in (|l8|) and (|2^) leaves the direction of the field 
vectors unchanged. Consequently, the result that they are orthogonal to each other and to the ray does not depend on 
the metric. But in the general metric the light ray is not a geodesic and consequently the vectors parallel transported 
along the ray will point in another direction than the physical ones. Thus what makes the metric (|40| ) special is the 
fact that the parallel transported and the physical field vectors coincide along the whole ray. 

Note that the energy densities, like the Poynting vector defined later, are time-averaged quantities with the fast 
temporal variation being disregarded. It should be mentioned that the numerical values of the energy densities (jsj) 
and (^5|) depend on the metric chosen. This is due to the fact that the physical quantity involved is the integral of 
the density over some volume, but as the line element is transformed according to 

cfe" = idx" (56) 



one has to change the density as well in order to keep the integral constant. Eq. (53) suggests that we introduce the 
auxiliary fields 

A^cr = y/eeo£cr and Ma ^ ^fi^ioHa- (57) 
These vectors have equal absolute value and fulfill the simplified equation 

2V'^/Co.V/35 = -ICaV'^VpS. (58) 



Eqs. ( [46| ) and (58) do not only describe the propagation of field vectors along the ray but they are also the key to 
understanding the propagation of the field intensity. The intensity / of a light ray is defined as the absolute value of 
the Poynting vector 

= ^v"''''Spn^ (59) 
and it is easily shown that it can be written in the form 

I^2^We. (60) 



Inserting (^6|) into (|60|) and taking into account that e/i = 1 we get for the intensity transport the result 

y^'SVpI^^.f^V^SVpieSj'^) (61) 



= i . l^{2e£''V'^SVi3£a + EaV^SV pe) (62) 
2 V Mo 

= -\.[^[e£'' £a{V pVf^S + V^SV p\ne + V'^SV plrip)] (63) 
2 V Mo 

= -IVpV'^S. (64) 

The right-hand-side of (Q) describes the change of the electromagnetic intensity along the ray. Due to the fact that 
/i and e are inverse, the terms containing these quantities cancel and the final result does not depend on their actual 
values. Consequently, light propagation in the geometrical approximation does not crucially depend on whether Ji 
and £ are constant or not, and neither the trajectory nor the energy flow are dependent on that. Clearly, the change 
of the fields E°' and H" along the ray does explicitly depend on the variation of /2 and e, respectively. The purpose 
of the last term in eq. (|4q) is, in a sense, to adjust the field in a way that, despite the changing permeabilities, its 
intensity flow is half of (^J). As has been pointed out before the product of Ji and e can be seen as an effective index 
of refraction. Consequently, the mathematical framework of the changed metric leads to a constant velocity of light 
and therefore corresponds to an empty but possibly curved space, for all physical processes where only the index of 
refraction is involved. 

The result in eq. (|6^ ) describing the change in intensity along the ray sheds light on the wave structure of the 
discussed set up. This can be easily understood if one looks not just at one trajectory but rather considers a thin 
but finite bundle of trajectories. Clearly, the distance between the trajectories may depend on the curve parameter. 
In the homogeneous case, however, all light rays are straight lines, and consequently the dependence is linear making 
the second derivative of the distance vanish. One can easily show that in this case the Laplacian of the phase d"daS 
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is proportional to the mean curvature, the arithmetic mean of the two principal curvatures, of the wave front. The 
formula for the change of the element of surface area spanned by a bundle of straight lines is given in books on 
elementary differential geometry (see e.g. [|5|) where the term "parallel surfaces" is used to denote the system of 
successive wave fronts. Consequently, the change in intensity is governed by the curvature of the wave fronts, whereas 
the value of the intensity is inversely proportional to the area of the cross section of the bundle. In the inhomogeneous 
case the right-hand side of eq. (^8|) cannot be interpreted in this simple way however because one has to add a term 
containing a Christoffel symbol to the simple Laplace expression. The additional term can be traced to the fact that 
in the inhomogeneous medium geodesies do not continue straight on, but rather curve through the dependence on 
the structure of the medium - in general in a different way at different points. Usually, the distance between these 
geodesies does not depend linearly on the curve parameter, i.e., the second derivative of the distance does not vanish. 
This effect is known as the geodesic deviation and is a characteristic of a curved space. Consequently, the changed 
law of intensity variation can be seen as a sign of the fact that the inhomogeneity of the medium creates an effectively 
curved space. 

More sophisticated mathematical methods should allow for a discussion of curved surfaces in a curved space. In 
this way (64) might turn out to be the curved-space analogue of the mean curvature. 



V. THE RIEMANN-SILBERSTEIN VECTOR IN INHOMOGENEOUS MEDIA 



By analogy with the fields Ma and JCa in cq. ( p7| ) one can generally introduce the fields 



seo 



(65) 



and 



(66) 



These fields can be seen to be the real and imaginary parts of the Riemann-Silberstein vectors 



(67) 



They allow for an alternative description of electrodynamics with the complete Maxwell equations in vacuum taking 
the form 



z— f - ±cV X F± 
ot 



and 



V • F± = 0. 



(68) 



Clearly, only one of the fields is needed for a complete description of electrodynamics. When used to describe 
photons, the Riemann-Silberstein vector has to be a correctly defined wave function, i.e., it may only include positive 
frequencies. In this case both vectors, {F+)a and (F_)q,, have to be used, corresponding to light which is left- and 
right-circular polarized, respectively, i.e., they then correspond to positive and negative helicity states In an 
inhomogeneous medium the field equations for the Riemann-Silberstein field become considerably more complicated 
than eq. (|68|). Written in the general form with an arbitrary diagonal metric they attain the form 







' dt 








' dt 





i((^,),V,lnyf + (F_),V,ln-L= 



1 



^\{F^)^^f,\nJ^ + {F+)^Wp\n 



1 



= 
= 



(69) 



and 



VaF_^ = ^ fev„ln^ 
2 V ^ ^ 



F^V„ln 



VaF^ = i (i^^V^ln^ 
2 V V^M 



(70) 
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Note that the two helicity states usually couple, i.e., helicity is not a constant of motion due to a term proportional to 
the derivative of In \fiije- In the non-physical case discussed before, i.e., when e and /x are proportional, the equations 
obviously decouple and one gets for the wave equation in the transformed metric the same form as in (|8|): 

-^ ^'g^2^" +i?.^^^ + V^Vff(F,)^=0. (71) 

As has been emphasized before, this choice of parameters allows us to describe light in an inhomogeneous medium 
exactly as if it were propagating in a curved space. In this case all polarization states are well defined and conserved. 
When as the ratio of e and [i becomes variable, the equations (69) and ( |70| ) couple again. But still, in the geometrical 



approximation both fields propagate independently, because the wave equation takes the form: 

- V^V^F_, + VpF-^V'ln — - V,F2V^ln — - 0. (72) 



Here all terms have been neglected which in the geometrical approximation contribute to the highest order in 1/fco 
only, i.e., terms not containing any derivative of F". Note that in this approximation the two helicity states decouple 
in the wave equation in all metrics. In the case where e and /i are inverse the two last terms on the right-hand side of 
([T^) vanish, leaving us with a standard wave equation expressed in terms of covariant derivatives. The discussion of 



geometrical optics in this case is then fully equivalent to that given in Section [V . Consequently, if the metric (|4C|) is 
applied, the complex vector J^" - defined analogously to - is parallel transported along the ray, which is a geodesic 
of the introduced metric. This means that the new metric allows to introduce well-defined polarization states which 
are propagated along the ray. In particular, helicity is a conserved quantity. Thus, if the Riemann-Silberstein vector 
had to be used to describe the wave function of a photon, the extension of this concept to inhomogeneous media 
should be discussed in the framework of quantum optics in a curved space. 



VI. SUMMARY 



In the present paper a general geometrical description of electromagnetic phenomena in inhomogeneous dielectric 
media has been given with the emphasis being on light propagation in the approximation of geometrical optics. To 
find the conceptually simplest description we introduced a variation of the spatial metric allowing for all metric tensors 
connected to the unit matrix by a conformal transformation. The equations for the fields and for the propagation of 
light rays can easily be rewritten in this general geometry when care is taken to transform correctly the fields and 
the material parameters. It turns out that when the metric tensor is equal to the unit matrix multiplied by - the 
square of the index of refraction - the equations become particularly simple: light rays are geodesies of the metric and 
the field vectors are parallel transported along the ray. Even the wave equation simplifies significantly. As long as fj, 
and £ in the standard metric are independent quantities the system in the transformed metric corresponds to a curved 
space filled with a dielectric medium, but with the transformed quantities /2 and e having a constant product. This 
property leads to significant simplifications in the description of light propagation. In the special case where /i and e 
differ only by a constant factor the situation becomes even simpler with the transformed qualities becoming constants. 
In this - non-physical - case the transformation fully eliminates the medium and all electromagnetic effects appear as 
if they were taking place in an empty but curved space. Whereas the wave equations contain significantly more terms 
when an effective medium is present, there are only slight differences between the two cases in the approximation of 
geometrical optics. The results there depend mostly on the behavior of the effective index of refraction which becomes 
constant. 

The results presented here contribute to a better fundamental understanding of electromagnetism in dielectric 
media and show the power of geometrical concepts in classical fields of physics. In particular the mathematical 
techniques of the general theory of relativity turn out to be well-suited for calculations in inhomogeneous media. In a 
forthcoming paper the presented results will be extended to electromagnetic fields in moving inhomogeneous media, 
thus supplementing the calculations presented in our papers on light in moving media. A tempting problem might 
be a geometric theory of anisotropic media which will require more sophisticated mathematical methods than those 
presented here including the introduction of a more general space than the Riemannian space. Another mathematical 
challange would be a discussion of the wave fronts in the spirit of differential geometry of curved surfaces in a curved 
three-dimensional space. Furthermore our results show that the correctly defined Riemann-Silberstein vector keeps 
a well-defined helicity during its propagation through the medium. This might trigger further investigations in the 
search for a photon wavefunction in a dielectric medium. 
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